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We report the thermodynamic properties of an ideal boson gas confined in an infinite periodic 
array of channels modeled by two, mutually perpendicular, Kronig-Penney delta-potentials. The 
particle's motion is hindered in the x-y directions, allowing tunneling of particles through the walls, 
while no confinement along the z direction is considered. It is shown that there exists a finite Bose- 
Einstein condensation (BEC) critical temperature Tc that decreases monotonically from the 3D ideal 
boson gas (IBG) value Tq as the strength of confinement Po is increased while keeping the channel's 
cross section, a^ay constant. In contrast, Tc is a non-monotonic function of the cross-section area 
for fixed Po- In addition to the BEC cusp, the specific heat exhibits a set of maxima and minima. 
The minimum located at the highest temperature is a clear signal of the confinement efi'ect which 
occurs when the boson wavelength is twice the cross-section side size. This confinement is amplified 
when the wall strength is increased until a dimensional crossover from 3D to ID is produced. Some 
of these features in the specific heat obtained from this simple model can be related, qualitatively, 
to at least two different experimental situations: '*He adsorbed within the interstitial channels of a 
bundle of carbon nanotubes and superconductor-multistrand-wires NbaSn. 

PACS numbers: 03.75.Hh; 05.30.Jp; 05.70.Ce; 67.40.Kh 



I. INTRODUCTION 

The properties of quantum systems at low dimension- 
ality and temperature have attracted the attention of 
researchers for a long time. In particular, the study 
of phase transitions, such as superfluidity, superconduc- 
tivity, or Bose-Einstein condensation, has been marked 
by an impressive interest among scientists in the field 
even though "true" long-range order phases are excluded 
by the Mermin-Wagner-Hohenberg theorentii^ in dimen- 
sions lower or equal than two. 

In two dimensions, however, the superfluid or super- 
conductor transition arises by the acquisition of quasi- 
long-range-order in the system as described by the 
Kosterlitz-Thouless phase transition^il. This, has been 
experimentally confirmed in thin superconducting—"— 
and helium-filmsii"— . Historically, different experimen- 
tal methods have been set up to study low dimensional 
phase transitions since the discovery of helium super- 
fiuidity. Indeed, studies on the behavior of helium 
films have been performed by using several experimen- 
tal techniques, like the adsorption of helium on simple 
plane substratesi^, or even on more complex ones as in 
nanoporous media such as cylindrical pores of Anoporei^, 
Vycor— or GelsiU^ glasses, where the effects of the sub- 
strate structure on the specific heat has been reported 
among other properties. More recently, studies of adsorp- 
tion of atoms and molecules on planar substrate o^^i^^ ; the 
two-dimensional character of high critical temperature 
superconductivity; and the discovery of graphene, have 
made two-dimensional systems to be widely explored con- 
trary to the case of quasi-one-dimensional ones. 

Although the theoretical aspects of one-dimensional 
systems have been extensively studied^S, only very re- 



cently ID experimental reports have attracted a great 
attention. For example, now days it is possible to cre- 
ate a ID Bose gas in cigar-shaped magneto-optic traps^i 
where the particle density, the cigar size and the inten- 
sity of the interaction between particles are experimen- 
tally tunable parameters, allowing one to examine quan- 
tum phenomena such as the superfluid to Mott-insulator 
phase transition o^^i^^ . 

On the other hand, with the advent of carbon nan- 
otubes, the realization of phase transitions in quasi- 
one-dimensional systems of different substances adsorbed 
on nanotube bundles is now possible^ir— . The quasi- 
one-dimensional character of these structures is a con- 
sequence of the enormous aspect ratios that nanotubes 
exhibit, with cross-sections in the nanoscale regime. On 
such length scales the single-particle energy levels cor- 
responding to the cross-section degrees of freedom are 
"frozen" leading to effective one-dimensional systems. 

In some other recent theoretical studies^i^, the oc- 
currence of BEC of a weakly-interacting quantum gas of 
Bose particles (parahydrogen or ^He) adsorbed within 
the interstitial channels (IC) of a bundle of poly-disperse 
carbon nanotubes has been predicted. The reported BEC 
transition and particularly, the dependence of the spe- 
cific heat on temperature, exhibit features of four dimen- 
sions in contrast to the expected one-dimensional behav- 
ior that has been, indeed, observed in the experimental 
report of Lasjaunias et al^ of the specific heat of ad- 
sorbed ^He in nanotubes. The authors of the former ref- 
erences, Refs. j28ll29| . argue that the presence of nonuni- 
formity in the nanotube cross-section gives rise to three 
additional degrees of freedom (the radii of the three tubes 
that form the IC), needed in their analysis to ensure the 
occurrence of the BEC transition which, as they claim^i 
based on the results reported in Ref. [s^] , doesn't exist in 
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a uniform bundle of identical one-dimensional nanotubes. 

Truly, it is well known that there is no BEC of a non- 
interacting boson gas in an impenetrable onc-dimcnsional 
box potential. Also, that there is no BEC whenever 
the spatial dimensions of at least one direction is finite, 
thus excluding the possibility of BEC in just one chan- 
nel. Therefore, a collection of independent IC's can not 
develop a BEC unless a coupling mechanism, such as 
tunnel-effect, is present between adjacent channels. A 
possible mechanism that couple the different IC's is con- 
sidered in Refs. f28ll29l |. where the authors argue that 
such a coupling leads to an effective density of states 
from which a non-vanishing BEC critical temperature 
is obtained. This effective density of states can be de- 
scribed as an inhomogeneously broadened convolution of 
the density of the heterogeneous transverse states with 
the one-dimensional density of states of a free particle 
that moves along the nanotube axis. 

Here we suggest that the packing defects due to the 
heterogeneity of the nanotube cross-section^ may lead to 
a tunnel-cffect-like mechanism that couples the different 
interstitial channels in the bundle, thus making BEC pos- 
sible. Indeed, it is well known that suitable confinement 
potentials can make the BEC long-range-order character 
to be stable against long-rangc-fluctuations. Such is the 
case of the ideal Bose gas in two dimensions trapped by 
harmonic potentials, which undoubtedly undergoes BEC. 
So, even in the case of a collection of nanotubes with uni- 
form cross-section, one would expect a BEC transition 
at a finite critical temperature when tunneling between 
neighboring channels (modeled using penetrable poten- 
tial barriers) is considered. 

In addition to the theoretical interest, the model 
presented in this paper can be applied to estimate 
critical temperatures and thermodynamic properties 
of superconductor-multistrand-wircs whose technologi- 
cal applications, that go from Nuclear Magnetic Reso- 
nances to magnets used in high-energy accelerators, im- 
pel their understanding. For example, several authors 
report experimental studies on the thermodynamic prop- 
erties of multistrand NbaSn (an Af 5 type of superconduc- 
tor) wires, that exhibit a transition temperature around 
18 K and can support fields up to f5 Tesla^. Typical 
numbers of filaments range from 10^ to 10^ NbaSn su- 
perconductor wires with diameters varying from a few to 
tens ^m. On the other hand, NbaSn multifilament wires 
come in different sizes, shapes and compositions, depend- 
ing on the techniques used to create them. The most 
common arc the Bronze Route, the Internal Sn diffusion 
process and the Powder Metallurgy (PM) methods, and 
they may have a core either of Cu, Sn, NbSn or NbCu 
alloys, while usually immersed in Cu, and the use of tan- 
talum and/or titanium barriers to prevent the Sn pollu- 
tion in Cii^. Experimentalists report either total spe- 
cific heat curves and/or curves that subtract the normal 
state specific heat, to avoid phonon and unpaired elec- 
tron interferenc e'^ '^''^^1'^^ . The main feature in this curves 
is the transition around 18 K with a typical width of 5 



K. However, in some cases a second peak appears around 
9 K which the authors interpret as the trasition of the 
remnants of unreacted Nb. Here, in light of our results 
we suggest an alternative interpretation to the meaning 
of these maxima. 

In this paper we report the thermodynamic properties 
with emphasis in the BEC critical temperature and the 
specific heat, of an ideal Bose gas within an infinite pe- 
riodic array of tubes. Our results are benchmarks for 
ongoing studies on the properties of real spatially con- 
fined systems such as: a) He atoms in interstitial carbon 
nanotube bundles, b) Cooper pairs in periodic tubes like 
multistrand NbaSn bundles or Bechgaard-salts^ or c) 
bosonic atoms in two dimensional opto-magnetic traps^^. 

Although at very low (or zero) temperatures and/or 
high densities the interaction between particles cannot 
be neglected, we focus on an interactionlcss boson gas 
to study the effects of a periodic confining potential on 
the properties of the system. We show that this sim- 
ple model captures qualitatively the properties of real 
systems, including the emergence of thermal phase tran- 
sitions and/or dimensional crossovers^^i^. 

In the following section we describe our system model. 
In Sec. HI we calculate the Bose-Einstein condensation 
critical temperature in addition to the specific heat and 
other relevant thermodynamic variables. In Sec. IV we 
discuss the results and present our conclusions. 



II. PERIODIC TUBE BUNDLES 

Our system model consists of N non-interacting bosons 
confined in an infinite periodical array of penetrable 
tubes of rectangular cross section of sides Qx and Qy, and 
infinite length. We model the tubes array by consider- 
ing two perpendicular Kronig-Penney (KP) delta barriers 
in the x and y directions with no constraints in the re- 
maining z direction (see Fig. [T]). Our periodic structure 
resembles either the bundle of homogeneous nanotubes, 
the superconductor-multistrand-wires or the experimen- 
tal 2D periodic lattice of tightly confined potential tubes 
created in Ref. [1^ . If the interaction between bosons is 
ignored, the Schrodinger equation for each boson of mass 
m in this system is 

Vix,y)} il}{x,y,z) = ekip{x,y,z) (1) 
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where Vx and Vy are the delta strength in the x and y 
directions, respectively. 

The Schrodinger equation ([1]) is separable in each di- 
rection such that Ek = Sk^ + £ky + £k^ is the energy per 
particle, where 
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with kz = 2'Knz/L the momentum in the z-direction, 
~ 0, ±1,±2,... due to the periodic boundary condi- 
tions in a box of length L, and Sk^ ^ are imphcitly ob- 
tained from the equations^ 




m[(aoiai)~i sm ap^a, + Pi {sin aoiflj -cos aoiai)(aoiai)~^] 
with aoitti = \j2mepilW and vahd for > 0.06. As ex- 
pected, Mi grows when Pi increases and for Pi > 40 the 
relation is almost linear. Even though ([5]) might seem 
to be a good approximation to calculate the thermody- 
namic properties of the systent^ in the low temperature 
regime, there are distinct effects, particularly in the spe- 
cific heat, that can only be observed when the full band 
spectrum is considered as is shown below. 




FIG. 1: (Color online) Periodic array of square cross section 
tubes. 



{Pi/aiat)sm{aiai) + cos{a,ai) = cos(feiai), (4) 
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rewrite the dimensionless constants Pi = mviai/h as 
P, = {mviXo /h'^){ai/Xo) = Poi(a,/Ao), where Ao = 
/i/VSTrmfesTo is the de Broglie thermal wavelength of an 
ideal boson gas in an infinite box at the critical temper- 
ature To ^ 27rn2„^/Vw'^BC(3/2)2/3 ~ ^Mfi^n^j^^ /mke 
with ub = N / the boson number density and the 
distance between the deltas barriers along the i = x and 
y directions. Poi is a measure of the tube wall impenetra- 
bility directly related to the delta-barrier strength. This 
can be understood if we recall the transmission coeffi- 
cient for an one dimensional delta potential of strength 
Vi, Ti = 1/(1 -|- Pq^/E)^ with Poi as defined above and 
E is the energy, in t? j^ma^ units, of particles arriving 
at right angles with the delta potential. We recover the 
following two limits: when Ppi goes to infinity the trans- 
mission coefficient vanishes and our model becomes an 
infinite number of decoupled tubes; when Pq^ = 0, = 1 
and the confining tube walls disappear recovering the 3D 
ideal Bose gas. By performing a series expansion of the 
left-hand-side of ^ just above of the exact Pi-dependent 
single-particle-ground-state energy eoi in the first band, 
the single-particle energy spectrum can be written to first 
order as 



III. CRITICAL TEMPERATURE AND 
SPECIFIC HEAT 

A. Grand potential 

The thermodynamic properties of the system are ob- 
tained from the grand potential il{T, , fi) for a boson 
gas 

n{T, L^,n) -TS - fiN 
= r!o + fcBr^ln[I-e-^("'=.-^)] (6) 

where ^Iq is the contribution of the ground state ki ~ 
with i = X, y, z, and /3 = l/kgT. As usual, U, S and fi 
denote the internal energy, the entropy and the chemical 
potential respectively. By using the dispersion relations 
given by ([3]) and ([4]) and after some algebra we obtain 



n(T,L^,n) = fcijrin[l 



L^m^'^ 1 



dk, &y53/2(e-^(^''^+^'^«-'^)), (7) 



where we have replaced the summations by integrals 
— >{L/2Tr f J d^k, assuming h^/mL'^ < ksT , 
and we have introduced the Bose functions^ ga{t) = 
X^z^i^V^'^- So = £ox + £oy is the ground state energy 
which depends on Poi and on Ui/Xo. 

From ([7| the thermodynamic properties for a 
monoatomic gas can be calculated using the relations 
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where z = exp(/3/^) is the fugacity. 



£0i 



■(1 — cos kio). 



(5) 



This is the dispersion that has been considered by sev- 
eral authors^i^ in the standard nearest-neighbor hoping 
approximation of the well-known Hubbard model with 
Epi = 0. The effective mass Mi is given by the expression 



B. Critical temperature 

We define the critical temperature Tc as the tempera- 
ture when the number of bosons in the ground-state level 
ceases to be negligible, i.e., Nq{Tc) ~ and the chemical 
potential fi{Tc) fio = fJ-ox + fJ-Oy = eo- 
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From the first expression in Eq. 
obtain the particle number 



and Eq. we 



N = 



1 
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At T = Tc, the first term vanishes and the critical 
temperature is obtained from 
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(10) 



Here we set the number density N/L^ equal to that of an 
ideal Bose gas (IBG) in the thermodynamic limit, with a 
BEG critical temperature Tq. Note that all the integrals 
involving the energy-spectrum in the x and y directions 
can be split in a sum of integrals over the energy bands 
folded in the first Brillouin zone. 

In the isotropic case, where Pqx ~ Poy = Po and 
o-x = OLy ^ flj the critical temperature as a function of 
the parameter Pq is shown in Fig. [2] for different values 
of the tube cross-section. Note that as the impermeabil- 
ity Pq of the tube walls increases, the critical temperature 
diminishes monotonically from Pq. In contrast, the vari- 
ation of Pc as function of a/Ao shows a non-monotonic 
behavior (Fig. [3]), decreasing for a/Ao < 1 and increasing 
otherwise . 
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FIG. 2: (Color online) Critical temperature in units of Po as 
a function of Po for different values of a^/Xo ~ tty/Ao = a/Ao. 

For both Figs. [2] and [H we find a similar qualitative 
behavior of the critical temperature as that reported for a 
boson gas in multilayers^^^, namely a decrease in Pc/Pq 
as Pq increases and a trend of Tc/Tq to go back to unity 
as a/Ao increases after having reached a minimum. How- 
ever, we notice that for similar Pq and a values as those 
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FIG. 3: (Color online) Critical temperature in units of Po as 
a function of ax/Xo = Uy/Xo = ffl/Ao, for different values of 
Po. 



used in multilayers, we obtain even lower Pc/Po values 
for bosons in square cross-section tube bundles , showing 
that the presence of an additional KP delta potential in 
the two-dimensional array emphasizes even more the ef- 
fects of confinement. Although we have only shown the 
critical temperatures for the isotropic ease, we will show 
in the following sections, the effects of anisotropy in the 
internal energy and specific heat by considering a rect- 
angular cross-section. 



C. Internal Energy 



by 



The temperature-dependent internal energy U is given 
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(e-£o)5i/2(x) + ^ff3/2(x) , (11) 



where e = e^^ -f- Sky and x = e.^^^^^^^ ■ 

In Fig. 0] we show the internal energy referred to the 
system ground state energy, for a square cross section, 
o-x = fly: a-nd in Fig. [5] we show results when Ox is differ- 
ent from Oy. In the isotropic case, the internal energy 
curves increase monotonically for Ox^y larger than Aq. 
However, when the separation is ax,y/Xo — 0.5 or lower, 
the effect of the KP potentials is revealed by two maxima 
that become more pronounced as the tube cross-section 
decreases. This behavior of the internal energy is simi- 
lar, albeit more complex, in the anisotropic case where 
the maxima and minima of the corresponding isotropic 
eases are still present. Minima are associated with the 
loss of freedom of particles in two directions, which we 
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FIG. 4: (Color online) Isotropic case. Internal energy as a 
function of T /To, for tube arrays of square cross section of 
several side sizes a^/Xo = a-y/Xo and Po = 100. 



D. Specific heat 

From Eq. ([8]) and (fTTj) . the specific heat becomes 

Cv I mP 



NkB N V &1T^h? Jo Jo 



dK dky [gi/2ix) 
+2/3(e - eo),9-i/2(x)(e - M + T^) + ^.93/2(x) 



(12) 



For T < Tc the chemical potential = /io is a constant, 
d^i/dT = and, using xo = e^'^^^^^''^ the last equation 
for the specific heat becomes 
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+I3{e - Mo)5-i/2(xo)(e - Mo) + ^53/2(xo)| • (13) 
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FIG. 5: (Color online) Anisotropic case. Internal energy as 
a function of T/To, a^/\o = 0.5, Oy/Ao = 0.1 and Po = 100 
(black line.) For comparison, the isotropic curves are also 
shown (color lines). 



will discuss in more detail when we analyze the specific 
heat behavior. 
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FIG. 6: (Color online) Isotropic case. Specific heat in Nks 
units, as a function of T/To, for different Ux = ay values and 
Po = 100. 

In Fig. [6] we show the specific heat for the isotropic 
case, as a function of T/Tq for different values of = 
Qy = a and Po fixed to 100. The BEC critical temper- 
ature of the system is displayed as a sharp peak in the 
specific heat that occurs at temperature whose value is 
lower than the critical temperature of a free ideal Bose 
gas. Except for the case with a^/ Xq ~ 0.1 where two min- 
ima appear, the specific heat shows only one minimum at 
a characteristic temperature that reveals the only length 
scale of the system, namely ax ^ ay. Such temperature 
satisfies the relation A ~ 2ax and marks the maximum ef- 

1/2 

fects of confinement, where A = (2nh^ /mkBT) is the 
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FIG. 7: (Color online) Specific heat in Nks units, as a func- 
tion of T/T(), for P() = 5000, a^^/Xo = ffly/Ao = 0.5 and 0.1 
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a, as is shown in Fig. [71 In both cases, Qx/Xq = 
dy/Xo — 0.5, and 0.1 the specific heat Cv/Nks, ap- 
proaches the one-dimensional classical value 1/2 over a 
relatively large region of temperatures. This behavior is 
more pronounced as Pq is increased. 

In Fig. [5] we can sec that for the anisotropic case 
where Ox o,y, the BEC transition occurs at a tempera- 
ture between the respective critical temperatures for the 
isotropic cases a = max{a^,ay} and a ~ min{a^,aj,}, 
and that the minimum corresponding to each isotropic 
case, both appear in the anisotropic one. 
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FIG. 9: (Color online) Specific heat in Nks units, as a 
function of T/Tq, for a^/Ao ~ 0.5 and Uy/Xo = 0.1 and 
Po — 100 (black line) compared to their respective isotropic 
cases: ax/Xo = Uy/Xo = 0.1 and 0.5, (color lines). One band 
only. 



E. Density of states 



FIG. 8: (Color online) Anisotropic case. Specific heat in Nks 
units, as a function of T/Tq, for a^/Xo = 0.5 and Uy/Xo — 
0.1 and Po ~ 100 (black line) compared to their respective 
isotropic cases: a^/Xo = ay/X^ = 0.1 and 0.5, (color lines). 



boson thermal-wavelength. There is another characteris- 
tic temperature for which the boson thermal wavelength 
is less than or equal to approx. Q.7 Ox that marks the 
point at which the effects of the walls are less conspicu- 
ous, above this temperature the system shows the usual 
3D free IBG behavior. The onset of this 3D behavior is 
observed as the last maximum on the right in the specific 
heat curve. 

As Pq becomes larger, the specific heat of the isotropic 
case reveals a one- dimensional behavior in a particular 
range of temperatures determined by the length scale 



The effects of the band structure are conspicuously ex- 
hibited in the density of states (DOS) 

.9(e) = X! ^ ^{^-^k.- Sky - EfeJ, (14) 
which can be written in the thermodynamic limit as 

^^'^ " (2^ E / , / , 

/oo 
dk^ (5(e - £fc^ - Efc^ - efcj, (15) 
-oo 

where we have explicitly written the integrals over the 
energy-spectrum in the x and y directions as a sum over 
bands of the integrals over kx and ky in the first Brillouin 
zone. 
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Upon integration over dk^ we obtain 
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where is the Heaviside step function. For energies 
close to the minimum, g{t) varies as e^/^ as it does the 
DOS of a free particle in three dimensions. This can be 
shown by noting that for energies close to eoi and for 
small fc, expression ^ can be approximated by + 
lT?kf/2Mi. In the isotropic case [Exq = e^o = eo/2 and 
Mx ~ My ~ M ), we can therefore write 

L3 



9{e) 



(27r)3 
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where e = e — £o, fc^ 



fcj + fc,^ and k is the maximum 

J, y 

value of fc in the first Brillouin zone. The exact value 
is not needed as long as e < [h?/2Ma'^) fc^. Thus, after 
evaluating the integral by a change of variable we obtain 



2to\^/^ 2M 



(17) 



(27r)2 \ ) h? 
The DOS g{e) from Eq. ^ is plotted in Fig. [10] for 
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FIG. 10: Density of states for tubes with Po = 100, a/Ao = 0.1 
(full dots) and a/Ao = 0.5 (open dots). 

Po = 100, ax/Xo = tty/Xo = O.I and 0.5. For ener- 
gies around the bottom of the first band given by p7|) is 
shown in the inset where a smaller value of Pq has been 
used to reveal the 3D behavior. The Pq dependence of 
g{e) enters only through the effective mass M, and in 
general, g{e) is a rather complex function of the energy 
(full-doted curve), however, in the limit of large values of 



Po the reminiscent behavior e~^^'^ of the one-dimensional 
DOS of a free particle is observed as a "falling" (see for 
instance the curve in open-circles) of g{e) with e. This 
behavior occurs when the energy-states along the z di- 
rection are the only ones that contribute to g{e), the con- 
tribution due to the other two directions being a constant 
related to the energy band-gaps. 



IV. DISCUSSION AND CONCLUSIONS 

At this point let us make the following remarks. In 
Fig. [9] we plot Cv as a function of T/Tq for different 
plane separations and Pq = 100, where only the first en- 
ergy band, exactly computed from Eq. ^ in the x- and 
?/-directions, has been considered. This situation qualita- 
tively resembles the results obtained if the dispersion re- 
lation ([5]) were used^i^. We consider the corresponding 
calculation of Cy in Fig. [51 where we have used enough 
number of exact energy bands needed to achieve the re- 
quired precision in the calculation of this property. We 
note that both curves coincide quantitatively and qual- 
itatively for temperatures not greater than Tq; for tem- 
peratures between To and about 5 To the qualitatively 
behavior remains best for the smaller plane separations; 
finally, for temperatures greater than 5 Tq even the qual- 
itative agreement disappears. Important structural in- 
formation is missing when just one band is included in 
the calculation, for instance, the second minimum that 
appears in the isotropic case of Fig. [8| (full-black line) 
is missing in the corresponding curve in Fig. [3) Clearly, 
the use of the lowest band leads to a one-dimensional 
behavior in the classical limit even though the system is 
three-dimensional. 

In the isotropic case (that can be related to homo- 
geneous nanotubes bundle-system), our model predicts 
a one-dimensional behavior of the specific heat in a 
range of temperatures determined by the distance be- 
tween the delta-barriers. Such results shown in Fig. [7| 
can be used to comparatively explain the one dimen- 
sional character of the specific heat of adsorbed *He in 
single-wall nanotubes as expected and observed in ex- 
perimental situations'^. Although Lasjaunias et al. did 
not report the presence of the peak that marks the BEG 
phase-transition we explain this fact in terms of the finite 
size of the experimental system where such transition can 
be made much less conspicuous. 

One important result is that our system always ex- 
hibits a EEC which is possible due to the coupling be- 
tween the different channels the bosons move in. When 
the wall impenetrability goes to infinity our system be- 
comes a set of decoupled tubes with zero EEC critical 
temperature as expected. Some results, often found in 
the literature, which claim that EEC in a bundle of ho- 
mogeneous nanotubes is not possible^, agree with ours 
if homogeneous means non-communication among tubes. 
However in the experimental set up, it is clear that the 
heterogeneity of the interstitial channels may have strong 
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effects on the thermodynamic behavior of the system, in 
particular in the BEC critical temperature which can be 
different from zero if the heterogeneity of the channels 
causes exchange of atoms among them. 

On the other hand, the use of multifilamentary su- 
perconducting tapes^ or wirea^ has improved the coil 
performance to support higher critical current density 
useful to create higher magnetic field to be used, for ex- 
ample, in the Large Hadron Collider currently under op- 
eration or the International Thermonuclear Experimen- 
tal Fusion Reactor planed to work in 2019. In both cases 
tapes and wires gather many filaments where pairs flow 
preferentially along the longitudinal direction. Critical 
temperature distribution in the NbaSn strands as well 
as the specific heat have been reported. The authors 
point out two transitions in the specific heat curve: one 
at the critical temperature (around 18 K) they associate 
to the complex superconductor wire and another one at 
a lower temperature (around 9 K) they associate to the 
unreacted Nb. However, if the Cooper pairs should be 
considered as bosons, our specific heat calculations show 
at least two critical temperatures: the lower one associ- 
ated to the collective effect of the filaments (our tubes) 
and a smooth maximum corresponding to the boson gas 
behavior in the individual filaments. In other words, the 
meaning given to the maxima in our Fig. [8] and in the 
Fig. 4 of Ref. [S^, are interchanged. A way to elucidate 
this controversy would be to make wires with larger fila- 
ment diameters, then prove that their smooth maximum 
shifts to the left as it is observed in our calculations. 

We summarize our main results in the following list: we 
observe that in the presence of periodical structures con- 
structed with orthogonal Dirac combs inside an infinite 
box filled with bosons, the critical temperature decreases 
from the 3D ideal boson gas Tq as Pq increases, while the 
plane separations ax,y/Xo are kept constant. It becomes 
zero when the periodic delta potential strengths become 
infinite. In other words, there is not BEC critical tem- 



perature different from zero for bosons in a tube of finite 
cross section and infinite length. 

As the separation between planes is lowered, the criti- 
cal temperature reaches a Po-dependent minimum value 
and then it is expected^i^ to increase again towards Tq. 

For systems with a^^y > Aq the numerical calculations 
for the critical temperature and specific heat arc very 
sensitive to the number of energy bands considered. To 
attain convergence we need to include up to 1000 bands. 

At T = Tc, the specific heat is continuous but has a 
discontinuity in its derivative. In the isotropic case it has 
one minimum and one or two maxima. The minimum is 
associated to particle trapping between two planes when 
its thermal wavelength is equal to 2a. This is corrobo- 
rated in the anisotropic case where the specific heat shows 
not one but two minima associated with the particle trap- 
ping in the x or y directions. 

The maximum at higher temperatures is associated to 
the onset of the system's approach to a 3D IBG behavior 
in this regime where the thermal wavelength A < 0.7a. 

While there is still a controversy over whether or 
not Bose-Einstein condensation of ^He exists inside in- 
terstitial filaments in bundles of carbon nanotubes, we 
conclude that in order to have BEC there must be a 
way through which the interstitial channels are coupled 
among them, either by effects of inhomogeneity of the 
tube bundles or by tunneling across the weaker intersti- 
tial walls. 

Finally, we mention that the proposed model in this 
paper gives account of systems composed of a very large 
number of quasi-one dimensional systems such as: bun- 
dles of carbon nanotubes, superconductor- multistrand 
wires, Bechgaard salt or 2D opto-magnetic traps. 

We acknowledge the partial support from grant PA- 
PIIT IN105011, IN117010-3, and IN106908. CONACyT 
104917. 
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